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ABSTRACT 

The CoRoT field LRa02 has been observed with the Berlin Exoplanet Search Telescope II 
(BEST II) during the southern summer 2007/2008. A first analysis of stellar variability led to 
the publication of 345 newly discovered variable stars. Now, a deeper analysis of this data set 
was used to optimize the variability search procedure. Several methods and parameters have 
been tested in order to improve the selection process compared to the widely used J index for 
variability ranking. This paper describes an empirical approach to treat systematic trends in 
photometric data based upon the analysis of variance statistics that can significantly decrease 
the rate of false detections. 

Finally, the process of reanalysis and method improvement has virtually doubled the number 
of variable stars compared to the first analysis by Kabath ct al. A supplementary catalog of 
272 previously unknown periodic variables plus 52 stars with suspected variability is presented. 
Improved ephemerides are given for 19 known variables in the field. In addition, the BEST II 
results are compared with CoRoT data and its automatic variability classification. 

Subject headings: binaries: eclipsing — methods: data analysis — stars: variables: general 
Online-only material: color figures, figure set, machine-readable and VO tables 



1. INTRODUCTION 

During the last decade, ground- and space- 
based surveys have been very successful in de- 
tecting transiting exoplanets. In addition to their 
primary science goal, the large photometric data 
sets acquired by them allow studying millions of 
stars for variability. Numerous projects thus pro- 
vide an exceedingly increasing number of detec- 



tions that are collected by variable star catalogs 
such as the General Catalo gue of Variable Stars 
(GCVS; ISamus et al.ll2009l) or the Variable Star 
IndejQ (VSX). Such catalogs not only broaden the 
statistical sample of variable stars, but are also im- 
portant to gain further knowledge about the dif- 
ferent processes that cause stellar variability. 

Several methods have been proposed to search 
for periodic signals in astr onomical time series (for 
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a goo d overview, see, e.g.. ISchwarzenberg-Czerny 
19991 ). One of the most widely applied algo- 



rith ms is the analysis of varian ce (AoV) statis- 



tic ( Schwarzenberg-Czernv 19961) . which provides 



an optimal period search in uneven sampled ob- 
servations. It has b een used very suc cessfully 
by pr ojects like HAT dBakos et al.l l2004h. WASP 



e.g.. iMac icicwski et al. 1 120111 ). or OGLE (e.g., 



ISoszvnski et al. I 20081) . In addition, the J index 



([Stetson' 1996 ) is frequently used to quantify vari- 
ability in general and/or for selecting candidate 
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; (e.g., IZhang et al. 


2003: 


Pasternacki et al. 


2011 


)• 



However, both methods - the AoV period 
search and Stetson's variability index - arc 
strongly affected by systematic trends present in 
ground-based data s e ts (se e , e.g.JPepper fc Burke I 



2006; Karoff et al. 



2007: Kabath et all 12009s 



Hartman et al. 1 120111) . Most dominant are diur- 
nal systematics, introducing artificial variability 
with periods of one day or multiples thereof. Such 
trends generally yield a higher ranking of non- 
variable stars, thus increasing the false alarm rate. 
A common approach to account for candidates 
with systematic variability is to set limits, e.g., to 
exclude detections within certain period ranges. 
However, any such manual mechanism is usually 
not well applicable to other data sets or projects, 
and the number of missed detections (false neg- 
atives) is often unknown. Therefore, a more so- 
phisticated treatment of systematic variability in 
combination with period search and ranking is 
needed. 

The Berlin Exop lanet Search Teles cope (BEST; 



Rauer et al. 120041) and BEST II ( Rauer et al 



20101 ) are used to perform groun d-based sup- 



port of the CoRoT space mission (jBaglin et al 
200G). By obtaining high-precision and long-time 
series photometry of the CoRoT target fields prior 
to the satellite's observations, planetary candi- 
dates can qu ickly be checked in the BEST data 
archive (e.g.. iDeee et al. 20091 : [Rauer et al.ll201dt 



Alonso et al.l l2012t ). In addition, the obtained 
light curves can be used to identify variable stars. 
The BEST project has already yielded the de- 
tection of several hundreds of new periodic vari- 
able s tars (iKaroff et al.ll2007tlKabath et al.l l2007. 
l2008l l2009bl lal: IPasternacki et al. 11201 ll ). 

A first characterization of periodic stellar vari- 
ability in the CoRoT fie l d LRa0 2 has been pub- 
lished by iKabath et all (|2009al) . Following the 
detections from this first publication (denoted as 
Paper I hereafter), it was possible to optimize 
the BEST II reduction pipeline. The large, well- 
characterized data set LRa02 was used as a proxy 
to analyze and automatically include systematic 
biases in the variable star candidate selection and 
period search, yielding a significantly reduced false 
alarm rate. 

The optimized procedure to rank periodic vari- 
able stars is described in this paper. In addition, 
we present a large catalog extension to the vari- 
able star classification in LRa02 that was obtained 
through the reanalysis and compare our results to 
publicly available CoRoT data. 



2. TELESCOPE AND OBSERVATIONS 

BEST II is located at the Observatorio Cerro 
Armazones, Chile. It is operated by the Institut 
fur Planetenforschung of the Deutsches Zentrum 
fur Luft- und Raumfahrt in robotic mode from 
Berlin since summer 2007. 

The system consists of a 25 cm Baker-Ritchey- 
Chretien telescope with a focal ratio of //5.0, 
yielding a wide field of view (FOV) of 1?7 x 1?7. 
It is equipped with a 4k x 4k, 16 bit Finger Lakes 
Imager CCD (KAF-16801E1) with a pixel size of 
9 /Ltm, and an angular resolution of l'/5 pixel -1 . 
BEST II observes without any filter to maximize 
the photon yield - the CCD sensitivity peaks at 
650 nm and is roughly comparable to the Johnson 
R band. 

The CoRoT long-run field LRa02 was observed 
by BEST II for 41 nights from 2007 November 
to 2008 February prior to the satellite observa- 
tions. As the FOV of CoRoT is slightly larger 
than BEST II, we split the field into two sub- 
fields (called LRa02a and LRa02b, respectively) 
and pointed at them alternating. Paper I indi- 
cates their corresponding center coordinates and 
shows the orientation with respect to the CoRoT 
FOV (Figure 1 in Paper I). 

3. DATA SET AND VARIABILITY CRI- 
TERIA 

The acquired data set was calibrated and re- 
duced using the BEST aut omated phot o metri c 



(|2010h . 



data pipeline as outlined in iRauer et al. 
The actual reduction procedure for the field LRa02 
has already been described in detail in Paper I and 
below only those steps relevant to the new results 
are discussed. 

In Paper I, Kabath et al. reported the detec- 
tion of 350 periodic variable stars (of which five 
were previously known) . Thei r selection was based 
upon the variability index J ( Stetsonl|l99 6), cal- 
culated for each star by 



En— 1 
k=l w k 



(1) 



where k is indexing individual data points. Pk 
is calculated from each pair of subsequent magni- 
tudes rrik and TOfe+i using the corresponding nor- 
malized residuals 4 and 4+1 ^ 



Pk = with 4 = 



n f mk — m 



n-l V (T fe 



(2) 

where <ik denotes the uncertainty of measurement 
k, m the mean magnitude, and n the number of 
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measurements for the selected star. The weights 
Wk in Equat ion Q were calcul ated following the 
approach o f IZhang et all (|20Q3l) as 



w k = exp 



At 



(3) 



where is the time of observation k and At is the 
median of all pair time spans (tfe+i —tk)- 

In Paper I, the limit of J > 0.5 was ap- 
plied in order to distinguish variable from con- 
stant field stars. This preselection yielded 1,858 
stars in LRa02a and 1,868 stars in LRa02b, re- 
spectively. For each star, th e AoV statistic 9 
( Schwarzenberg-Czernv 19961 ) was calculated for 
a period range of 0.1-35 days. The light curves 
were folded with the frequency w m ax correspond- 
ing to the maximum AoV value, 



B(w max ) = max(6(w)) , 



(4) 



and then inspected visually. Most folded light 
curves showed no clear periodic variability or an 
artificial period of one day or multiples thereof, 
which can be caused by systematic effects due to 
the observational cycle. iKabath et al. ( 2009al ) fi- 
nally identified 173 periodic variables in LRa02a 
and 177 in LRa02b. 

4. REANALYSIS 

In the BEST archive, some stars are marked as 
clear variables with large J indices (up to 10 and 
higher) . Figure [T] shows the count of all stars in 
subfield LRa02b and the corresponding number of 
variable star detections in Paper I as a function 
of the J index. Altogether, the large number of 
false positives shows that the J index alone is not 
an effective criterion for selecting variable stars. 
Furthermore, a number of clear detections with 
low J values indicates that several variables must 
have been missed in Paper I due to the cutoff. 

In order to improve the quality of our preselec- 
tion process and to maximize the detection yield, 
we performed a deeper analysis of our data set 
LRa02. Field LRa02 was chosen for this pur- 
pose because it was observed with longest duration 
within the BEST project so far. 

The performed reanalysis consisted of three 
main steps. 

First, the best period was searched in all light 
curves from the initial data set without any prese- 
lection criterion, i.e., directly using the AoV mul- 
tiharmonic algorithm for all stars in both subsets. 
Each of the resulting 37,361 folded light curves in 
LRa02a and 66,974 in LRa02b was then examined 
visually for periodic stellar variability. In addition 



to the variable stars from Paper I, this first step 
already revealed 189 additional variable stars that 
were not detected in our first publication. Dur- 
ing this step, we also discovered a bug in our im- 
plementation of the J index that yielded system- 
atically lower values especially for short periods, 
which is why many of these detections show peri- 
ods of less than a day. 

Second, the combined data set of variable stars 
from Paper I and the additional manual detections 
were used to optimize the BEST II selection pro- 
cess (see next Section[5]). With a very good knowl- 
edge of this data set, it was possible to compare 
different selection methods and to adjust their cor- 
responding parameters. 

Third, the new search algorithm was applied 
to the data set with optimized parameters. The 
improvements to the pipeline finally lead to an in- 
crease of detections by another 135 variable stars 
that were not found during all previous steps (see 
Section O). 



5. IMPROVEMENTS ON VARIABIL- 
ITY SEARCH 

The deep analysis of the BEST II data set 
LRa02 gives us the opportunity to study the per- 
formance of variability search algorithms in detail. 
The aim is to recover all variable stars in the data 
set automatically and to minimize the number of 
false alarms (and manpower) at the same time. 

In this section, we show the limitations of the 
J index with regard to systematic trends (Sec- 
tion I5.ip . describe how to quantify the perfor- 
mance of a variability search fSection !5.2p . present 
the algorithms tested (Sections 15.31 and I5.4|) . and 
finally compare the performance of different ap- 
proaches and parameters (Section 15. 5p . 

5.1. Limitations of the Variability Index J 

After the first step of the reanalysis, the vi- 
sual inspection, the majority of new detections 
showed J indices below the limit of 0.5 applied 
before, which is why they were not detected in Pa- 
per I. The relation between periods and variability 
is shown in Figure [2] A clear bulk of stars can be 
found at about 0.01 < J < 0.1, with the limits be- 
ing widely period independent and populated by 
stars showing noise only. Most new variable star 
detections are found in the region between J = 0.1 
and the cutoff limit of J — 0.5 from Paper I. A 
small number of new variables with J > 0.5 were 
not detected in Paper I because their J indices 
have been initially underestimated (see Section[4j. 

The dominant variation in many light curves is 
due to diurnal systematics, aliasing, or a combi- 
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Fig. 1. — J index histogram for subfield LRa02b (qualitatively equal to LRa02a). The solid line shows the total star count in 
bins of 0.1 (J < 2) and 1 ( J > 2), whereas the dotted line shows the number of variable star detections from Paper I. The red 
dashed line denotes the cutoff limit of J = 0.5 as applied in Paper I. 




0.1 1.0 10.0 

Folded period [d] 

Fig. 2. — Variability J index plotted vs. determined pe- 
riod without treatment of systematic effects for all stars 
in the BEST II field LRa02a (top) and LRa02b (bottom). 
Variable stars identified in Paper I are marked with blue 
dots, whereas variable stars from the manual reanalysis 
appear red. The dashed line shows the selection limit of 
J = 0.5 as applied in Paper I. 



nation of both. Figure [2] shows a large accumula- 
tion of stars having periods of one day or integral 
fractions/multiples thereof, often in combination 
with large J indices. Consequently, this leads to a 
very high number of false alarms when using the 
Stetson index as the only criterion for variability 
selection. In the example of the data set LRa02, 
a cutoff limit of J = 0.1 would be sufficient to 
include all variable stars in the selection sample, 
but only 74% of all light curves would be sorted 
out. The remaining large sample of 31,000 stars is 
mainly affected by systematic effects and contains 
only 681 stars with real physical variability (see 
Section |6|) . The corresponding false alarm rate of 
about 98% shows the need for an automated treat- 
ment of systematic variability, which is not part of 
the J index. 

5.2. Quantitative Assessment of Period 
Search Algorithms 

Two fundamental criteria are used to assess the 
quality of period search algorithms: the signifi- 
cance of the detection itself and the correct deter- 
mination of the frequency of variability. 

First, a quantity £ is introduced to evaluate the 
detection efficiency of any given search algorithm. 
Detection methods are usually based upon a single 
numerical value q (e.g., the Stetson index, q = J) 
that can be used to prioritize a candidate list. The 
success of ranking variable stars high in the list is 
measured with £ for each tested search algorithm. 
It ranges from for the perfect algorithm (all pre- 
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Fig. 3. — Schematic view on the rejection of systematic 
frequencies. The set £1* (711,712) of non-systematic frequen- 
cies is obtained by excluding the subsets S7^ ya (master power 
spectrum cut) and f2| ys (empty phases). The size of both 
can be adjusted with the parameters ni and ri2- Maxima of 
0(oj) and SQ(ui) are searched within Q*(ni,n2) to obtain 



the frequencies it>m,L and 



,(2) 



respectively. 



viously identified variable stars listed first) to 1 
(listed last). For details on the calculation of £ see 
the Appendix. 

The second criterion is tested by comparing the 
frequency w C orrcct that was verified manually with 
the frequency of a tested algorithm. We consider 
a tolerance range of 2% around w correc t for a cor- 
rect determination. Also included are 2% devia- 
tion around half or twice that value, because the 
distinction between these is often ambiguous from 
the light curve itself. The fraction n u of correctly 
identified frequencies can then be used for a quan- 
titative comparison between tested algorithms. 

5.3. Frequency Determination and Exclu- 
sion of Systematics 

Because systematics and their aliases are usu- 
ally limited to a set of few well-defined frequen- 
cies {w sys }, they can be excluded by searching 
the best frequency CL> max only on a subset tt* = 
{uj*} = {u>}\{uj a y S } (Figured). We tested three 
different methods to account for systematic fre- 
quencies, both independent of each other as well 
as in combination. 

1. Master power spectrum. Systematic periodic 
signals affect many light curves in a data set and 
can thus be distinguished from real stellar vari- 
ability by analyzing many power spectra 0j(w) of 
individual stars i statistically. We use the mean 
of all spectra to build a master spectrum 



(5) 



i=l...N, 



ter spectrum, a baseline fit is determined. Among 
several investigated functions, a polynomial log- 
log fit, 



in(e$rM) 



i=0 



(lnw)*, 



(6) 



with coefficients Cj maps the baseline best and 
most reliable for various tested data sets. The 
baseline-subtracted spectrum can then be searched 
for systematics. Figure [4] shows how Om — 0^f sc 
peaks clearly at the diurnal frequencies that we 
aim to identify. Finally, a simple cutoff using the 
standard deviation ct^| so of the subtracted spec- 
trum 9m — 6^f sc is applied to filter systematics 
automatically. The corresponding set of frequen- 
cies having peaks at least m • o-\j sc above the 
average power spectrum is defined by 



base l 
T M h 

(7) 

where the parameter n\ can be adjusted to quan- 
tify the degree of exclusion. 

2. Frequencies with empty phases. Ground- 
based observations are strongly affected by pe- 
riodic gaps in the data, most commonly due to 
the diurnal cycle. The incomplete phase coverage 
leads to aliasing and can often cause false posi- 
tive detections. We use a simple model to exclude 
frequencies with poor phase coverage: first, the 
folded light curve is split into N boxcs = 100 in- 
tervals of the same length. The number of empty 
intervals N^£J(lj) is then counted for each sam- 
pled frequency. A frequency is considered system- 
atic by this criterion if the fraction of empty inter- 
vals is larger than a relative threshold parameter 
7i 2 , i.e., 



nT(n 2 ) = {uj\NZZ t s y (u J )>n 2 -N hc 



.}■ (8) 



Both criteria [7] and [8] are merged to exclude 
systematic frequencies determined in either way, 
i.e., the overall set of non-systematic frequencies 
(Figure [3) is defined by 

o*(m,n 2 ) - M \ (nf s ( ni )nnT(n 2 )) (9) 



and the search for uj„ 



is restricted from the 



whole frequency range (Equation (|4)) to the sub- 
set Q*(ni, na): 



V W e^(n 1 ,n 2 ):9H<e( W « x ). 



(10) 



Note that the function 



N?™?*? (ui) is similar 



boxes 



but not equivale nt to the window function 7at(cl>) 
( Deeming"! 19751) . Limiting the latter was tested 



In order to identify significant peaks in the mas- 



as an alternative criterion; it shows a very simi- 
lar ranking performance, but slightly less correctly 
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Fig. 4. — Master power spectrum for data sets LRa02a (top) and LRa02b (bottom). 



determined frequencies. Being more simple, the 
empty phase criterion was chosen for the final test 
setup. 

3. Power spectrum scaling. In addition to the 
exclusion of systematic frequencies, we also inves- 
tigated a method to include the information about 
systematics into the AoV results directly. Instead 
of searching for the maximum of Q(ui), an artificial 
spectrum 



69(u) = e(u)/6jif(u) 



(11) 



is created by dividing the AoV spectrum of ev- 
ery star by the master power spectrum &m- Its 

(2) 

maximum is found at the frequency u>mix in anal- 
ogy to Equation (fT0|) : 

6 n*(ni,na) : 80 (u) < 50(^1) . (12) 

5.4. Variable Star Ranking 

In addition to the J index, two methods to pri- 
oritize variable star candidates were tested. 

The first method takes the AoV result directly, 
i.e., 

q<P = 6(^1). (13) 

In its special case of no excluded systematic fre- 
quencies (m — ► oo, ri2 = 1), this is a widespread 
method for prioritizing variable star candidates. 
Likewise, the maximum of the divided power spec- 
trum 50 could serve as a variability indicator: 



The AoV statistic 



e = (n - 



nu)\\xn || 2 



— n rra— COm- 

Jl|| ||X-X|||| 2 

pares the quadratic norm of a model xii (with 
nil free parameters) with the residuals that re- 
main after s ubtraction of the model from n ob- 
servations x (jSchwarzenberg-Czernv I [l999f h Be- 
cause it has optimum period d etection properties 
(jSchwarzenberg-Czernv 1996), we expect q% to 
yield the best ranking. However, as an empir- 
ical alternative we also tested the light curve's 
standard deviation with and without the periodic 
signal obtained by the AoV multiharmonic fit: 



a' (lu 



(i) \ 

max ' 



and 



(2) 



rr'(, , (2) ^1 
<7 [LO max ) 



(15) 



It is dependent upon the frequency w max deter- 
mined in the previous section, which is why differ- 
ent choices of u; max lead to different rankings. The 
quoted a' refers to the standard variation after 

subtraction of the co rresponding fit (for calcula- 

tion o f coefficients, see lSchwarzenberg-Czernv fc Kaluznv 
19981) . 



5.5. 



Comparison of Variability Search Per- 
formance 



The quantities £ and n, 
for both frequencies 



have been calculated 



„(i) 

^max 



and uj 



(2) 



and both 



(k) 

tested ranking methods q[ % . For each run, the 
parameters m and ni were varied independently 



qf = 8G(u,%&. 



(2) 



(14) 
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Table 1 

Results of the tested variable star ranking methods and parameters. 











(2) 

Wmax 










(2) 
01 




Index 




LRa02a 
LRa02b 


(0,0) (0,0) 

0.0051 0.010 

(0,0) (0,0) 

0.0035 0.0056 


(0.2, 10%) 

0.0035 

(0,0) 

0.0027 


(0.2, >70%) 

0.010 

(0,0) 

0.0062 


0.032 
0.034 




LRa02a 
LRa02b 


(5-10, 10%) 

65% 

(>10,0) 

63% 


(>5, 10%) 

84% 

(>5, 10%) 

81% 





Note. — The quantities £ and are shown for each tested ranking parameter q, both methods 

(k) 

of determining the best frequency Wmax and both analyzed data sets. For clarity, only the best 
value achievable by variation of the parameters m (master power spectrum cut) and ri2 (empty 
phases) is shown for each method. The corresponding values/ranges of n\ and 712 are shown in 
small brackets above each value. 



on the following values: 

m = {0,0.01,0.05,0.1,0.2,0.5,1,2,5,10,1000} 
n 2 = {0%, 10%, 20%, 30%, 40%, 50%, 
60%, 70%, 80%, 90%, 100%} 

The results of the comparison are summarized in 
Table [1] which shows the best value for £ and n u 
achievable with each tested method. 

5.5.1. Number of Harmonics 

Paper I and the first step of the reanalysis 
(Section |4j determined the stellar variability us- 
ing AoV periodograms with two harmonics, which 
was also used for this comparison. Furthermore, 
the number of harmonics was set to N — 7 in a 
second test in order to increase the sensitivity on 
sharp signals that are, e.g., caused by eclipsing bi- 
naries (for the sensitivity dependence on the num- 



ber o f harmonics, see, e.g.. lSchwarzenberg-Czernv 
1999). Both results show very similar ranking per- 
formances (for the best ranking method qi, we 
find £7 w £2 ± 0.001 in both data sets), but the 
N = 2 test naturally yields a slightly better (12%- 
15%) frequency match with the initial run that 
was obtained using the same number of harmonics. 
However, the test with seven harmonics revealed a 
number of additional interesting eclipsing binaries 
that could not be detected using the smaller num- 
ber of model parameters (see Sections 16 . 1 1 and HP)) . 
Therefore, we prefer the latter for our improved 
variability search (Section 15. 6p and focus in the 
following on the details of the search performance 
with N = 7 harmonics. 

5. 5. 2. Ranking 

The performance of the ranking differs only 
slightly between the tested methods. Figure [5] 



shows how the quantities qi (AoV) and q 2 (er- 
ratic) both provide a sorting that lists stars with 

real variability first. The numerical quantity £ 

(2) 

yields with 0.003 for q\ a minimum close to 
the optimal ranking (£ = 0). It is an order of 
magnitude lower than the previously used J in- 
dex (£ 0.03) and significantly lower than rank- 
ing the AoV power without exclusion of system- 
atics (£ 0.022 for q{ 1] with n x = 1000 and 
ri2 = 100%). In particular, this corresponds to 
a drastically decreased false alarm rate. 

Figure |5] shows the dependence of the ranking 
performance on the parameters n-y and n 2 for the 
two best methods q^ and q± . In both cases, 
the most restrictive exclusion of systematic fre- 
quencies yields the best sorting. Thereby, the 
cut in the master power spectrum (Equation ([7])) 
has a slightly larger impact than the exclusion of 
empty phases (Equation ©). The minimum of £ 
is reached for n\ = n 2 = 0, but is almost inde- 
pendent of ri2, because the first criterion is more 
restrictive. 

5.5.3. Frequency Determination 

The frequencies found in Paper I and by manual 
reanalysis of the data set are in good agreement 
with the values of Wmax (see Tableland Figure[7]). 
Without master spectrum division (fc = 1), about 
two-thirds of the frequencies are recovered. How- 
ever, the yield increases to about 80% if the pro- 
cedure is applied (fc = 2). 

Interestingly, the frequency exclusion from the 
first criterion (Equation ([7])) now has the reverse 
effect - the maximal agreement is reached if it is 
almost switched off by setting n\ > 5. Smaller 
values of ni are too restrictive and can increase the 
number of wrong periods by up to about 10%. On 
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Fig. 5. — Number of variable stars N' v as a function of the number of inspected stars N c * for LRa02a (left) and LRa02b 

(k) (k) 

(right). The different lines represent tested ranking methods: q\ , q^ , and the Stetson J index (left to right). Solid lines 
represent the unweighted case (k = 1), whereas dotted lines include the effect of master power spectrum division (k = 2). Only 
the parameters n\ and ri2 of the most successful sorting are used for each method (compare Table [TJ. The black dashed line 
shows the optimal ranking for comparison. 




Fig. 6. — Influence of systematic frequency exclusion method on variable star ranking efficiency for the data set LRa02b. The 
quantity £ is shown on the jy-axis as a function of the model parameters ni (lines) and (x-axis) for the best tested ranking 
method using q\ . The left plot shows the results with uj^I x as the maximum frequency, whereas the effect of division by the 
master power spectrum (Wmax) can be seen on the right. 




Fig. 7. — Fraction By of correctly identified frequencies uL„ as a function of model parameters ni (lines) and n,2 (a;-axis) for 
variable stars in data set LRa02a (left) and LRa02b (right), respectively. 
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the other hand, the exclusion of empty phases has 
again a small influence, although a value of n 2 = 
10% yields a slight improvement for the majority 
of tested scenarios (Figure [7]) . 

The remaining small group of variable stars de- 
tected with a different period has been analyzed 
carefully. The majority of them shows multi- 
period variation and was identified with a rational 
multiple (e.g., 1/7, 2/5) of the original frequency. 
For some stars, the original period had to be re- 
vised during the reanalysis (see also Section [6]) . 
A small rest shows amplitudes close to the noise 
level, such that the period could not be determined 
unambiguously. 

5.6. New Selection Method 

Based upon the results from the comparison, a 
new procedure was set up to search for variable 
stars within the BEST project. 

1. The J index is used to exclude non- variable 
stars in order to save computation time. 
High values of J can originate from either 
real variability or systematic trends, leading 
to a very high false alarm rate when being 
used as the only ranking criterion. How- 
ever, low J indices give a reliable criterion 
for non- variability, i.e., neither physical nor 
systematic variations. In the studied data 
set LRa02, no star shows clear variability 
below J = 0.1. This limit is therefore used 
for analyses of BEST II data sets, for which 
it typically excludes 50%-75% of all stars. 

2. The AoV algorithm is applied with N = 7 
harmonics to the selected subset (J > 0.1) 
in order to obtain power spectra. 

3. The improved selection method is applied to 
rank all investigated stars. Following the re- 
sults from Section [5. 5[ a master power spec- 
trum is calculated, the number of empty 
boxes is counted for each test period, and 
each individual power spectrum is divided by 
the master power spectrum (Equation (|1 1|) ) . 

(2) 

For the ranking, the frequency Wmax is de- 
termined from the subset of non-systematic 
frequencies in the divided spectrum SQ(uj) 
by following Equations ©-© and ([HJ with 
the parameters n\ = and ri2 = 10%. The 

( 2) 

corresponding maximum <50(wmax) (Equa- 
tion (|14|) ) serves as the quantity q for prior- 
ization (see example in Figure E]). However, 

(2) 

in order to improve the final period u)f, Wmax 
is recalculated without exclusion of system- 
atic frequencies from the master power spec- 



trum, i.e., by setting n\ — > oo, Ti2 — 10% and 
applying Equations j7])-(|9|) and (fT2|) again. 

4. All light curves are folded with their respec- 
tive final periods LOf and analyzed visually 

(2) 

in descending order of q = <58(cjmax)- 

6. RESULTS 

In addition to the 350 variables already pub- 
lished in Paper I, 279 stars in LRa02 were identi- 
fied with clear periodic variability (114 in LRa02a 
and 165 in LRa02b). Furthermore, we identified 
52 suspected periodic variable stars (17 in LRa02a 
and 35 in LRa02b). For the latter, the quality 
of the light curves is not sufficient to fully ex- 
clude systematic errors as sources of variability, 
the folded light curves are partly incomplete or the 
point-spread functions of two neighboring stars 
overlap each other such that the signal cannot be 
well separated. The number of detections in this 
work is compared with Paper I in Table [2j 

Table 2: Variable Star Detections in BEST II 
Data Set LRa02 - Sum marized Counts for Paper I 



( Kabath et al.ll2009d ) and This Work. 





LRa02a 


LRa02b 


Total 


Paper I 


173 (4) 


177 (1) 


350 (5) 


This paper 


114 (2) 


165 (5) 


279 (7) 


This paper (suspected) 


17 (0) 


35 (0) 


52 (0) 


Total 


304 (6) 


377 (6) 


681 (12) 



Note. The number of previously known variables in the 
field confirmed by BEST II is given in brackets (included 
in first number). 



6.1. Application of the New Method to 
LRa02 

A large number of 189 new variable stars was al- 
ready identified in the BEST II data set LRa02 by 
the first step of the reanalysis, the manual screen- 
ing (see Section |4j . After the search procedure 
was tested and optimized using the results from 
Paper I and this additional sample of detections, 
the most successful sorting method (see descrip- 
tion in Section [5]) was finally applied to search the 
data set LRa02 once more. 

In addition to the 350 variable stars published 
in Paper I and the 189 found by the manual re- 
analysis, this improved search yielded another 135 
previously unknown variable stars. Reasons why 
these went undetected by previous searches are as 
follows. 

• Systematic trends or aliases were found in- 
stead of the real periodicity. 
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LRo02b2_01822 



J=4.23 



LRo02b2_1 7835 



J = 0.24 




02468 10 02468 10 

u [days -1 ] u [days -1 ] 



Fig. 8. — Example for power spectra ©(a;) (upper plots) and <50(u) (lower plots, sec Equation {TTJ) - the star LRa02b2_01822 
(left plots) is strongly affected by systematic noise, while LRa02b2_17835 (right plots) shows both physical and systematic 
variability. In each plain AoV spectrum (upper plots), the position of the overall maximum is marked - it is found at 

systematic frequencies for both cases. Furthermore, the maximum of @(io) and the J index are both much larger for the first 
star, leading to a false alarm when using these quantities to rank the variability. The functionality of the new variable star search 
algorithm is shown in the lower plots. In addition to the division by the master power spectrum, systematic frequencies are 
filtered out (Equations (0— l(9j with parameters n\ = and 712 = 10%, marked red in the spectrum), and the maximum cjmix is 
determined on the non-systematic subset of frequencies (Equation J 1 2D ~) . To use the corresponding maximum = 50(cj^2 x ) 
for ranking is much more sensitive to real variability. The second star LRa02b2_17835 is a new detection of this paper (Table [3] 
and Figure [Tot . 




Fig. 9. — Photometric quality of the BEST II data set LRa02a (in line with LRa02b). The standard deviation a of each light 
curve is plotted vs. the BEST II instrumental magnitude Rg for each star. The red fit line indicates the lower noise limit 
o-g ln (iJg) for this data set. 
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• The AoV algorithm was run with N = 2 har- 
monics for Paper I and the manual reanal- 
ysis, but with N — 7 for the latest search. 
This leads to a higher sensitivity for non- 
sinusoidal variations, which is particularly 
important for the detection of Algol type 
eclipsing binaries (at least 27 additional EA 
variables with long periods (P > 2 days) can 
be attributed to the increase of N, some of 
them being very eccentric). 

• A total of 21 variable stars on the edge of 
the BEST II FOV with few datapoints were 
excluded by earlier reductions of the data 
set. 

Due to an improved data quality (Figure [9]) and 
increased sensitivity for non-sinusoidal events, we 
could also refine the periods for 17 of the variable 
stars published in Paper I (see Table 4). 

6.2. Classification of Variability 

The newly discovered periodic variable stars 
are classified by the shape, amplitude, and pe- 
riod of their folded light curves according to the 
scheme used by the GCVS. We identified pulsating 
stars of the Delta Scuti type (DSCT), RR Lyrae 
type (RR), and Cepheid variables (CEP). Eclips- 
ing binary star systems were classified as Algol 
type (EA), Beta Lyrae type (EB), or W Ursae Ma- 
joris type (EW). Stars having sinusoidal-like light 
curves and showing eclipses are marked as rotat- 
ing ellipsoidal variables (ELL), whereas some light 
curves exhibit characteristic features of spotted 
stars (SP). In case the type of variability could not 
be determined from the light curve alone, we clas- 
sified a star as VAR, meaning that further observa- 
tions are needed to constrain the proper physical 
origin of stellar variability. Stars that are variable 
on time scales longer than the observational base- 
line are classified as long periodic (LP). This class 
can also include non-periodic variables. 

6.3. Catalog of New Variables in LRa02 

The newly identified periodic variable stars of 
the BEST II data set LRa02 are listed in Ta- 
ble [3] A set of stars from Paper I with revised 
ephemerides is presented in Table |4] Due to the 
reanalysis, the internal numbering is not consis- 
tent with Paper I, which is why IDs were given 
different prefixes (i.e., LRa02a2 and LRa02b2, re- 
spectively) . Corresponding IDs of the nearest Two 
Micron All Sky Survey (2MASS) object are given 
if the catalog coordinates do not differ from our 
astrometry by more than 2 arcsec. In some cases, 
two neighboring stars could be spatially resolved, 
but their photometric apertures overlap so that 



the actual source of variation is ambiguous. Such 
objects have been marked with a "c" (crowded), 
whereas suspected variables are indicated with an 
"s" flag. All instrumental magnitudes given are 
obtained without any filter so that their absolute 
value should only be considered as a rough es- 
timationH Amplitudes and ephemerides (Tq in 
rHJD=HJD-2,454,400) are the results of the AoV 
algorithm. However, after visual inspection, sev- 
eral periods have been refined manually (usually 
to rational multiples of the initial result). 

The corresponding phase-folded light curves 
can be found in Figures \W\ and [TTJ For stars 
classified as LP, light curves are shown unfolded 
and no ephemerides are given in Tables [3] and 4. 
Light curves in a machine-readable format as well 
as finding charts are available upon request. 

6.4. Comparison with Known Variables 

We searched the VSX and the GCVS for pre- 
viously known variable stars within a radius of 10 
arcsec around each of our new detections. A total 
number of seven variables was found in these cat- 
alogs. For comparison, the new BEST II results 
for these are included in Table [3] and marked with 
"k" as known variables. 

The three stars CoRoT 110742676, NSVS 
12579155, and NSVS 12585233 have periods longer 
than 50 days. BEST II confirms their long-time 
periodicity, but the phase coverage of their cy- 
cles is insufficient to confirm the periods quan- 
titatively. For the four eclipsing binaries ASAS 
J064835-0534.3, DYMon, [KEE2007] 1318, and 
[KEE2007] 1334, both the classifications and pe- 
riods are completely confirmed. 

The latter two have first been detected by 



Kabath et all (|2007h in the CoRoT IR01 field with 
BEST, for which the chosen FOV shows a small 
overlap with the BEST II field LRa02a. However, 
the precision in the periods of these two binaries 
could be significantly improved because the LRa02 
data set covers a much larger time period (41 com- 
pared to 12 nights) and the photometric quality of 
BEST II is better. 

6.5. Comparison with CoRoT 

The LRa02 data sets of BEST II and CoRoT 
are, except for the fact that they point to the same 
field, completely independent of each other and 
thus provide a great opportunity to compare the 
scientific results of the two surveys. 



2 The zero point was shifted to minimize the residuals with 
the R magnitude of USNO-A2, yielding an rms deviation 
with the catalog of about 0.5 mag for all matched stars. 
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Table 3: Catalog of Variable Stars Detected in CoRoT Field LRa02 after Reanalysis of the Data Set, Sorted by Internal BEST II Identifiers. 



BEST ID Flag 


2MASS ID 


a(J2000.0) 


<5(J2000.0) 


Rb (mag) 


T [rHJDJ 


P (d) 


A (mag) 


J Index 


Type 


Other Names 


LRa02a 


LRa02a2.00759 


06473611-0352264 


06"47 r 


l 36.1 s 


-03°52'26' 


6 


18.11 


34.638 


0.828 ± 0.002 


0.4 ± 0.2 


0.306 


EA 




LRa02a2.01976 


06483833-0312255 


06 h 48" 


'38.3 s 


-03°12'25' 


6 


16.34 


35.163 


0.8348 ± 0.0006 


0.30 ± 0.05 


1.37 


EA/SP 




LRa02a2_03383 


06481122-0346301 


06 h 48 r 


'11.2 s 


-03°46'30' 


4 


17.13 


35.097 


1.294 ± 0.005 


0.12 ± 0.07 


0.323 


CEP 




LRa02a2.04701 


06472903-0432171 


06 h 47 r ' 


'29.0 s 


~04°32'17' 


2 


15.25 


38.772 


8.6 ± 0.2 


0.04 ± 0.02 


0.445 


ELL/SP 




LRa02a2.05225 


06483988-0339178 


06 h 48 r 


'39. 9 s 


-03°39'17' 


9 


16.51 


34.937 


0.4350 ± 0.0004 


0.12 ± 0.07 


0.812 


ELL/SP 




LRa02a2.05867 


06475987-0416581 


06 h 47" 


'59.9 s 


-04° 16' 58' 


3 


17.72 


34.841 


0.5255 ± 0.0005 


0.4 ± 0.2 


0.457 


EB 




LRa02a2.06108 


06475271-0424518 


06'' 47" 


'52.7 s 


-04°24'51' 


9 


17.00 


34.601 


2.747 ± 0.008 


0.22 ± 0.07 


0.469 


EA 




LRa02a2.06147 


06492768-0308050 


06 h 49 r 


»27.7« 


-03°08'05' 


2 


13.87 


34.711 


0.40998 ± 0.00006 


0.73 ± 0.03 


23.6 


EB 




LRa02a2.06168 


06483769-0348472 


06 h 48 r 


'37.7 s 


-03°48'47' 


3 


16.15 


35.256 


0.8230 ± 0.0006 


0.16 ± 0.03 


0.807 


EA 




LRa02a2.06363 s 


06482994-0356402 


06 h 48" 


'29.9 s 


-03°56'40' 


3 


15.80 


35.944 


1.646 ± 0.007 


0.04 ± 0.03 


0.548 


ELL 




LRa02a2.06993 


06482849-0403248 


06 h 48 r 


'28.5 s 


-04°03'24' 


8 


13.82 


34.705 


0.4248 ± 0.0002 


0.106 ± 0.008 


5.75 


EW/ELL 


CoRoT 300002950 


LRa02a2.07010 


06481669-0413092 


06 h 48 r 


"-16.7 s 


~04°13'09' 


3 


14.56 


35.923 


2.309 ± 0.006 


0.04 ± 0.02 


0.658 


EA 


CoRoT 300002493 


LRa02a2_07056 


06484660-0349153 


06 h 48 r 


'46.6 s 


-03°49'15' 


5 


17.20 


34.934 


0.5425 ± 0.0005 


0.22 ± 0.06 


0.838 


EW 




LRa02a2.07090 


06493379-0311214 


06 h 49 r 


'33.8 s 


-03°11'21' 


3 


17.24 


34.925 


0.5427 ± 0.0007 


0.15 ± 0.07 


0.374 


EW 




LRa02a2.07148 


06480066-0427255 


06 h 48 r 


'00.7 s 


-04°27'25' 


5 


16.74 


34.676 


0.35786 ± 0.00008 


0.60 ± 0.05 


1.62 


EW 





Notes. Stars have been matched with the closest 2MASS object within a radius of 2 arcsec. Coordinates are given for epoch J2000.0 and were derived by an astromctric match of CCD to 
USNO-A2 coordinates. Given magnitudes arc instrumental and reflect the CCD sensitivity as observations have been obtained without filter. Overlapping apertures of neighboring stars can lead 
to contaminated light curves — such cases are marked with a "c" flag. Suspected variables arc marked by "s". Amplitudes A and ephemcrides (variability period P and times of minimum 
brightness T 0j given in rHJD— HJD-2454400) are the results of the AoV algorithm. 

(This Table is available in its entirety in machine-readable and Virtual Observatory (VO) forms in the online journal. A portion is shown here for guidance regarding its form and content.) 



Table 4: Catalog of Periodic Variable Stars with Revised Parameters Compared to Paper I (jKabath et al.ll2009al) 



BEST ID 


Flag 


2MASS ID 


a(J2000.0) 


t)(./2000.0) 


Rb (mag) 


To [rHJDJ 


P(d) 


A (mag) 


J Index 


Type 


Other Names 


LRa02a2.04976 


k 


06473121-0432567 


06"47" 


'31 


2 s 


-04°32'56' 


9 


15.05 










0.964 


EA 


lra2a_00269 




LRa02a2.07847 


k 


06483188-0407577 


06 h 48" 


'31 


9 s 


-04°07'57' 


9 


16.45 


39.649 


1.5015 ± 0.0001 


1.3 ± 


1 


4.04 


EA 


lra2a_00416 




LRa02a2.08694 


k 


06491193-0342336 


06 fe 49" 


'11 


9 s 


-03°42'33' 


8 


16.14 


34.694 


0.31415 ± 0.00005 


0.46 ± 


03 


4.50 


EW 


lra2a.00450 




LRa02a2.10471 


k 


06500331-0315474 


oe^so" 


' 03 


3 s 


-03°15'47' 


4 


16.44 


35.159 


2.695 ± 0.001 


0.8 ± 


1 


0.943 


EA 


lra2a.00531 




LRa02a2.12019 


k 


06485506-0423379 


06 fe 48" 


'55 


I s 


-04°23'38' 


1 


15.33 










6.12 


CEP 


lra2a_00601, 


CoRoT 110826631 


LRa02a2.24433 


ck 


06520050-0326327 


06'' 52" 


' 00 


5 s 


-03°26'32' 


8 


13.66 


98.695 


6.785 ± 0.001 


0.6 ± 


1 


6.82 


EA 


lra2a_01126, 


CoRoT 110677259 


LRa02a2.24442 


ck 


06520077-0326255 


06'' 52" 


' 00 


8 s 


-03°26'25' 


7 


14.32 


98.695 


6.785 ± 0.001 


0.26 ± 


01 


1.69 


EA 


lra2a_01127 




LRa02b2.03306 


k 


06474385-0541125 


06 , '47 r ' 


'43 


8 s 


-05°41'12' 


7 


12.57 










7.77 


LP 


lra2b_01648 




LRa02b2.05167 


k 


06484535-0458261 


06 h 48" 


'45 


4 s 


-04°58'26' 


2 


15.69 


34.740 


0.7749 ± 0.0004 


0.43 ± 


04 


3.92 


EA 


lra2b_01600, 


CoRoT 300003628 


LRa02b2.13901 


k 


06484069-0537534 


06 fe 48" 


'40 


7 s 


-05°37'53' 


6 


14.33 


35.085 


0.52869 ± 0.00009 


0.39 ± 


02 


14.9 


EW 


lra2b_01437, 


CoRoT 300003441 


LRa02b2.23935 


k 


06491567-0550387 


06 h 49" 


'15 


7 s 


-05°50'38' 


9 


14.06 


35.241 


2.966 ± 0.009 


0.16 ± 


02 


2.53 


EA 


lra2b_01257, 


CoRoT 110657689 


LRa02b2.31945 


k 




06'' 50" 


'23 


9 s 


-05°26'29' 


5 


15.60 


34.697 


2.239 ± 0.004 


0.24 ± 


04 


1.42 


EA 


lra2b.01080 




LRa02b2_38262 


k 


06505364-0527405 


06 h 50" 


'53 


6 s 


-05°27'40' 


8 


12.43 


41.500 


41 ± 1 


0.18 ± 


01 


28.1 


SR 


lra2b_00968, 


CoRoT 110839568 


LRa02b2.48861 


k 


06512387-0549233 


oe^i" 


'23 


9 s 


-05°49'23' 


6 


13.58 










5.93 


LP 


lra2b.00738 




LRa02b2.51632 


k 


06511964-0605018 


oe^i" 


'19 


7 s 


-06°05'02' 


2 


15.92 


35.282 


1.388 ± 0.005 


0.08 ± 


05 


0.879 


ELL/SP 


lra2b_00687 




LRa02b2.59423 


k 


06522734-0545516 


06'' 52" 


'27 


3 s 


-05°45'51' 


8 


13.41 


34.730 


1.8214 ± 0.0001 


0.62 ± 


01 


14.4 


EA 


lra2b_00469 




LRa02b2.66851 


k 


06531961-0540496 


06 h 53" 


'19 


6 s 


-05°40'49' 


6 


15.29 


34.760 


0.45935 ± 0.00009 


0.57 ± 


03 


6.44 


EW 


lra2b_00323 






Fig. 10. — Phase-folded light curves of variable stars detected in field LRa02 after reanalysis of the data set. 
(The complete figure set (33f images) is available in the electronic edition of the Astronomical Journal) 



A step in the scientific analysis of CoRoT data 
consists o f an automatic stellar v a riabil ity clas- 
sification (jDebosscher et all 120071 l2009h . This 
method was also applied to the CoRoT observa- 
tions of field LRa02, and the results are mean- 
while - together with the full light curves - pub- 
licly available through the CoRoT archived 

The CoRoT LRa02 data set contains 11,448 
targets, from which 10,392 (91%) are matching a 
BEST II target within a maximum angular dis- 
tance of 1". Because LRa02 was not covered com- 
pletely by BEST II (Figure 1 in Paper I), 454 
CoRoT targets are located outside of the BEST II 
FOV. Furthermore, the magnitude ranges do not 
overlap completely, so that 425 bright CoRoT 
stars are saturated on the BEST II CCD. The re- 
maining 177 CoRoT targets are within the FOV 
and right magnitude range, but have no BEST II 
counterpart due to technical issues such as bloom- 
ing. In the same way, BEST II observed a total of 
93,943 stars in both pointings that have not been 
given a CoRoT mask or are located outside the 
CoRoT FOV. 

From the 681 variable stars presented together 
in Paper I and this work, 262 variables (190 from 
Paper I, 72 from this work) match a CoRoT target. 
The corresponding CoRoT IDs are shown for the 
new detections in Tables [3] and 4. 

CoRoT observed the LRa02 field about 1 year 
after BEST II, and the observing times do not 
overlap. Because this work focuses on the im- 



provement of variable star detection for ground- 
based telescopes, CoRoT and BEST II light curves 
have not been combined. However, we note 
that such a combination might yield improved 
ephemerides. BEST II light curves of interesting 
objects will be provided on request. 

6.5.1. Detection Efficiency of BEST II and the 
New Search Algorithm 

For stars both observed with CoRoT and 
BEST II, it is possible to investigate the per- 
formance of BEST II more in detail. The auto- 



3 CoRoT data are available to the community from the 
CoRoT archive: |http: / /idoc-corot. ias.u-psud.fr/ 



matic classification (jDebosscher et al. 1 120071 ) pro- 
vides information about the amplitude of vari- 
ation for CoRoT targets. Because the satellite 
has a much higher photometric precision, most 
of these amplitudes are well below the detection 
limit of BEST II. However, the knowledge of stel- 
lar variability with higher precision can be used 
to evaluate the detection efficiency of BEST II. 
Of particular interest are two questions: first, 
how many stars with sufficiently high variation in 
the CoRoT data set have been detected as vari- 
ables from BEST II data. Second, if these are 
clearly distinguished by the new detection algo- 
rithm from the rest of stars having variabilities 
below the threshold of BEST II. 

The full amplitude Ac of variation was derived 
from the fit coefficients given in the CoRoT clas- 
sification. As this contains each a low- and high- 
frequency entry for every CoRoT target, the am- 
plitude was calculated as the maximum of both. 
The value Ac itself contains no information on 
whether the signal can be detected by BEST II or 
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LRa02a2_04976 



LRa02a2_07847 



EA P = (T.5015 ± 0.0001)° 



LRa02a2_08694 



EW P=(0.3H15 ± 0.00005)° 




Fig. 11. — Phase-folded light curves of known variable stars with revised parameters. 
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Fig. 12. — Detection efficiency of BEST II for the field LRa02. Shown arc 10,392 stars that arc measured by both BEST II 
and CoRoT. The expected BEST II S/N of CoRoT amplitudes is plotted on the x-axis (Equation H16II ). On the y-axis, the 
left plot shows the suggested new quantity for variability ranking q = 8&(ui^l x ) (Equation 1141 ) and the right plot the Stetson 
J index for comparison. Variable stars detected in Paper I or this work are marked red. The area of possible BEST II detections 
is indicated in gray. 



not, which strongly depends on the magnitude of 
a given star. Therefore, the quantity 



S/N = A c /a'E iD (R B ) 



(16) 



is used to estimate the variability signal to noise 
(S/N). The noise a^ in (R B ) gives the photomet- 
ric precision achievable with BEST II in the given 
data set for a star of magnitude Rb- It was deter- 
mined by a fit to the a- magnitude plot of the field 
(Figure [9} . 

Figure [T2l shows the expected S/N for all stars 
that are contained in both data sets. From 680 
stars with CoRoT amplitudes that should be vis- 
ible in the BEST II data (S/N > 3), 162 were de- 
tected as variable stars in Paper I and this work. 
448 stars are expected to have a large S/N > 3, 
but show no significant variability in the BEST II 
data set (q < 9). 70 stars with both expected and 
measured variability (S/N > 3, q > 9) were not 
detected by BEST II. 

We have checked CoRoT light curves with sig- 
nificant variability in the CoRoT classification, 
but that are not variable in the BEST II data. 
Many of the m show strong instru mental effects 



(hot pixels; Auvergne et al. 2009) which obvi- 



ously mislead the automatic classification algo- 
rithm. However, from the 448 stars in this region, 
only 88 show a probability larger than 95% to be- 
long to any class so that most can be considered 
false alarms. 

Stars that are expected to be variable (S/N > 
3) and show variation in the BEST II light curve 
(q > 9), but were not detected as variable stars af- 
ter visual inspection have been re-inspected. Such 
targets have been missed due to the smaller phase 



coverage of BEST II, or because the real noise 
of individual light curves is underestimated using 
erg m (e.g., due to higher order extinction effects 
for very red stars). 

The distinction between variable and non- 
variable stars works very well using the new rank- 
ing quantity q = <59(wmax) (see Section [5]). If the 
limit q — 9 is chosen to separate variables from 
the bulk of non- variable stars, 92% of all matched 
stars are found below the limit. Only seven vari- 
ables have q < 9, but their light curves and low 
S/N indicate rather false detections than too low 
variability values. From all 780 stars with q > 9, 
one third belongs to the set of variable star de- 
tections. For comparison, Figure [12] also shows 
the J index vs. the expected S/N. The plot shows 
clearly that the separation between real and ar- 
tifical variability is much weaker. The strength of 
the new ranking is particularly clear in the regime 
of 1 < S/N < 3, i.e., close to the detection limit 
of BEST II. 

6.5.2. Comparison of Classifications 

The variable star detections of BEST II were 
compared in detail with the automatic CoRoT 
classification for the 262 matched variables. The 
overall agreement between both methods is very 
good; details regarding the period and magnitude 
determination as well as the classifications ob- 
tained by BEST II and CoRoT are given in the 
remainder of this section. Figure [13] shows some 
instructive examples of variable stars in agreement 
(a) and with differences in the determined periods 
and/or classifications ((b)-(f)). 

The mean magnitudes of matched stars are in 
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Fig. 13. — Examples for variable stars identified both by BEST II and CoRoT. Each row shows light curves of the same 
star: in the first column, the BEST II light curve is folded with the period from BEST II; in the second column, the CoRoT 
light curve folded with the BEST II period; and in the third column, the CoRoT light curve folded with the CoRoT period. 
Unfolded light curves are shown for long periodic variables. 
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Fig. 14. — Histogram of differences between CoRoT and 
BEST II magnitudes (for matched detections, without very 
long periodic variables). 
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Fig. 15. — Histogram of CoRoT and BEST II period ratio 
for matched variable star detections. The insets show the 
two main peaks enlarged (normalized). 



reasonably good agreement (Figure [14]) . Only a 
few very long-term variables - like the example 
in Figure [TSTd) - show differences in the order of 
1 mag and above because BEST II and CoRoT ob- 
served during different phases; the remaining ma- 
jority differs by only (0.059 ± 0.158) mag. 

Periods determined by BEST II have been com- 
pared with the main frequency from the auto- 
matic CoRoT characterization. For 72.5% of the 
matched stars, the periods are equal or integral 
(n = 1, . . . , 5) multiples of each other to a preci- 
sion of at least 1%. The histogram of period ratios 
in Figure [15] shows that most detections have been 
identified in the CoRoT data with half the period 
compared to BEST II. This is because many pe- 
riods are doubled during the visual inspection of 
BEST II light curves in order to show full cycles 
(e.g., Figures [TSTc) and ITBT e)). in particular for 
W Ursae Majoris eclipsing binaries (EW). 

The classes from visual inspection of BEST II 
light curves match the automatic classification of 
the CoRoT data set well. From the 262 variables 
present in both data sets, 196 stars have BEST II 
variability classes that are consistent with either 
the short- or long-periodic classification in the 
CoRoT data set. Note that the variability classes 
used by BEST II and CoRoT are slightly differ- 
ent: for example, BEST II distinguishes within 
the CoRoT class ECL between the eclipsing bi- 
nary types EA, EB, and EW, while the CoRoT 
scheme includes, e.g., slowly pulsating B-stars, 
which are simply identified as VAR within the 
BEST II study. All such consistent refinements 
are considered a classification match. 

The 64 stars with a clear disagreement in the 
variability classifications were checked carefully by 
reviewing the light curves from both BEST II and 
CoRoT. For 26 stars, the variability classes ob- 



tained by BEST II appear more realistic. Most 
of such cases are LP variable stars (e.g., Fig- 
ure E^d)) that have been identified as such in the 
BEST II data set by visual inspection. For these 
cases, even the longer CoRoT baseline does not 
cover a full cycle, so that the period and classifi- 
cation obtained from the CoRoT pipeline are not 
conclusive. Furthermore, some eclipsing binaries 
are clearly misclassified by the automatic CoRoT 
procedure. Figure 113( e) shows an example of a 
W Ursae Majoris type eclipsing binary that was 
identified as an RR Lyrae pulsator - most likely 
because it was detected with half of its physical pe- 
riod from the CoRoT data set. For three stars sim- 
ilar to the example in Figure fT3T b). narrow eclipse 
events were not detected by the CoRoT analysis. 
Most of the stars with implausible variability types 
were classified as BE by the automatic classifi- 
cation , which was described by iDebosscher et al. 
(2009) as a "trash" class regarding its wide pa- 
rameter spread. For 16 stars like the example in 
Figure [RTT f). the CoRoT classifications are in bet- 
ter agreement with the measurements. For most 
of these cases, this clearly results from the bet- 
ter photometric quality of the satellite data. For 
22 cases, the photometric data itself are insuffi- 
cient to choose between the CoRoT and BEST II 
classifications (e.g., EW/ELL). No automatic clas- 
sification data were available for the two CoRoT 
targets 110833621 and 300001413. 

7. SUMMARY AND DISCUSSION 

The CoRoT target field LRa02 was observed 
with BEST II during 41 nights from 2007 Novem- 
ber to 2008 February (see Paper I). We reanalyzed 
the data set in order to improve the detection 
method and to maximize the number of detec- 
tions. In addition to the 350 periodic variable stars 
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already published earlier (Paper I), we present a 
catalog of 272 new variables and 52 stars with sus- 
pected variability. For seven known variables both 
the classifications and periods are confirmed, and 
the periods could significantly be improved for two 
of them. Revised ephemerides are presented for 17 
variable stars from Paper I. 

From a manual inspection of all light curves 
without any preselection criterion, it turned out 
that most of the new periodic variable stars went 
undetected in Paper I because the applied vari- 
ability criterion using the Stetson J index was too 
restrictive. But because the J index is heavily bi- 
ased by systematic effects, a smaller cutoff limit 
leads to very high false alarm rates (w 98%) and 
is therefore not a practical alternative. However, 
although the J index is not capable to distinguish 
between systematic and stellar variability, it can 
still be used to exclude non-variable stars from 
the analysis: no variable star is falsely rejected if 
light curves with low variability indices J < 0.1 are 
sorted out in both data sets of this study. This sep- 
aration can be particularly useful if the full AoV 
process is too time consuming for a whole data set. 

The deep characterization of the data set en- 
abled us to compare and train different automatic 
methods for an improved variability ranking. In 
particular, a master power spectrum was calcu- 
lated as the mean of all individual AoV spectra. 
This method proved a valuable tool for exclusion 
of systematic frequencies and hence the ranking 
of real variability. The best algorithm found sepa- 
rates variable stars very effectively from the non- 
variable background population and in parallel re- 
covers their frequencies well. The new ranking 
method is particularly superior to the J index in 
regimes where the amplitude of variation becomes 
comparable to the noise level, i.e., close to the de- 
tection limit of the photometric system. It shows 
an almost equal performance for both independent 
subsets LRa02a and LRa02b, so that it should be 
easily applicable to other data sets. 

Finally, the results from both Paper I and this 
work have been compared with the publicly avail- 
able CoRoT data set. BEST II obtained light 
curves for 91% of all CoRoT targets in LRa02. 
From all 681 variable stars observed with BEST II 
in the field, 262 are matched with CoRoT targets. 
Stellar amplitudes measured by CoRoT show that 
BEST II detects variability efficiently in its pa- 
rameter range. The results of period determina- 
tion and variability classification are found to be 
in very good overall agreement, confirming the va- 
lidity of the measurements, data reduction, and 
scientific analysis of both projects. 196 stars show 
consistent classifications. CoRoT yields a better 



classification in case of small variability for 16 tar- 
gets, e.g., Be stars. On the other hand, the fi- 
nal visual screening of folded BEST II light curves 
proves to be an important step in the analysis. For 
26 stars it was found to yield a more realistic clas- 
sification, in particular for eclipsing binaries and 
stars with variability on time scales comparable to 
the observational baseline. 
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A. Calculation of sorting parameter 

Lot X* = {xi} with i <G {1, . . . , A*} be the group of all stars in a data set. The selection of variable stars 
- e.g., by sorting all stars by a ranking quantity q - can then be considered a permutation p v : X* — > X* 
that splits X* into a part X v C X* containing all A„ variable stars and another part with the rest: 

p v (xi) e X„ for i < N v 
p v (xi) g X v for i> N v . 

A check of the first N v within {p v {xi)} would thus reveal all variable stars in the data set. 

Unfortunately, such an optimal sorting p v is usually unknown. In practice, a given permutation p v aims 
at a similar splitting of variable and non- variable stars, but contaminates both groups with false positives. 
The number of identified variable stars N' v < N v thus depends on the number of stars N„ < that are 
actually checked: 



N' v { Pv ,N„) = Y^5 v {p v { Xi )) 

i=l 



with 8 v (xi) 



l Xi e x v 

otherwise 



A given variable star selection p v can be compared directly with the optimal procedure p v . The number 
of missed variable stars is 

N'v iPv, N„) = N' v {p v , N„) - N' v ( Pv , N„) , 

whereby 

N v (p v , N„) = mm(N„, N v ). 

For a comparison between different approaches, it is useful to evaluate the performance of p v as a whole. 
For that, one can define the quality parameter 

4 N V (N*-N V ) 

The parameter £ sums missed variable stars for all values of ./V c * . The denominator accounts for normalization 
such that £ = for p v (best selection) and £ = 1 for its counterpart (worst selection). 
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